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Abstract. A new class of solutions to the coupled, spherically symmetric 
Einstein-Maxwell equations for a compact material source is constructed. Some 
of these solutions can be made to satisfy a number of requirements for being 
physically relevant, including having a causal speed of sound. In the case of 
vanishing charge these solutions reduce to those found by Bayin and Tolman. 
Only the latter can be considered as having physically realistic properties. 
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1. Introduction 

One of the most fundamental problems in general relativity is the construction of 
exact solutions to the Einstein equations (or the Einstein equations coupled to other 
equations associated with fundamental physical theories). An even more challenging 
problem is is to find solutions that might represent physically realizable systems 
that could possibly occur in nature. It is well known that of the many solutions 
to the systems of equations derived from the Einstein equations lead to unphysical 
behaviour. In many cases, the spacetimes constructed may be truly singular even 
in the presence of matter or other fields. Alternatively negative energy densities 
and/or negative pressures are often required in order to meet regularity conditions or 
boundary conditions that are imposed upon the metric or other geometric objects. 
In other cases where energy densities and pressures are positive, often one finds that 
the speed of sound violates causality by being greater than the speed of light. For 
example constant density configurations, where the pressure has a radial dependence 
to ensure that the interior solution remains static will lead to infinite sound speeds. 
Charged dust solutions on the other hand have zero sound speed. They require such 
a fine balance between the mass and charge distributions that they are almost always 
unstable to any perturbations. They either undergo gravitational collapse or fiy apart 
due to electrostatic repulsion. 

This paper will present a new solution to the spherically symmetric time 
independent Einstein-Maxwell system of equations that govern the behaviour of the 
space-time in the interior of a charged fluid sphere. At the outer boundary the solutions 
will be matched to the external vacuum solution for the field equations i.e. the 
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Reissner-Nordstrom solution. It will be demonstrated that some of these solutions 
have parameter values that lead to physically realistic properties. 

While the total number of known exact interior solutions to the Einstein-Maxwell 
equations is much smaller than their uncharged counterparts, the number of solutions 
that can be said to be physically relevant is small. A review of such solutions has 
been provided by Ivanov [1] and this makes a good starting place for understanding 
the different approaches that have been taken and the results obtained from those 
approaches. In some instances corrections to some of the original papers arc provided. 
The most commonly used methods employ either isotropic coordinates (where the 
problem was first discussed by Papapetrou [2] and Majumdar [3]) or curvature (i.e. 
Schwarzschild-likc) coordinates. 

As mentioned in the introductory paragraph above, this paper will be primarily 
concerned with the construction of physically "interesting" solutions that might arise 
from regular initial conditions for charge and mass densities along with realistic 
pressure distributions. The motivations for this are to find compact charged fluid 
configurations that might represent realistic sources for the exterior vacuum Reissner- 
Nordstrom solutions. In particular one can ask whether it is possible to find realistic 
solutions for the extreme case where the total charge is equal to the total mass. 
A second motivation is to determine the possible end states for systems that have 
collapsed gravitationally but due to the repulsive nature of the electrostatic interaction 
are unable to form black holes and result in stable configurations. 

The approach to be taken here is to write the Einstein-Maxwell system in 
curvature coordinates. The Einstein equations will yield three equations for two metric 
functions in terms of the mass density, the charge aspect and the fluid pressure. The 
Maxwell equations simply state the the electrostatic Maxwell field is just the Coulomb 
field in a spherically symmetric space time. The solution will be given in terms of 
explicit closed- form functions of the radial coordinate for the two metric coefficients, 
a mass density, a charge aspect (the integrated charge density) and the pressure of 
the fluid. 

The review by Ivanov outlines that various ansatse that have been utilized to 
find exact solutions. Essentially one gives two out of the five functions and uses the 
Einstein-Maxwell equations to construct the remainder. Of course some methods of 
solution are easier than others, particularly if one is not concerned with ensuring 
that the solution have a physical interpretation. The simplest method is to give the 
functional dependence of the metric functions and construct the remaining variables 
from the derivatives on the "left-hand-side" of the differential equations. On the other 
hand giving the mass and/or charge distributions along with an equation of state of 
the form P = P{p) requires solving a set of nonlinear coupled differential equations. 
But the method has the advantage that there is some control over the physics. Other 
methods use a combination of the two approaches. 

The requirements that make a solution physically relevant has had much 
discussion in the literature. (See e.g. the reviews by Delgaty and Lake [4] and 
Finch and Skea [5] that analyze over one hundred known solutions for uncharged, 
spherically symmetric, fiuid solutions in general relativity.) Those that we expect to 
meet arc the following: 

(i) The metric coefficients are regular everywhere including at the origin r = 0. 

(ii) The metric functions match to the Reissner-Nordstrom functions at the fluid- 
vacuum interface. 
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(iii) The mass density is positive and decreasing outward toward the boundary. 

(iv) The integrated charge and mass increase outward to give the RN parameters at 

the boundary. 

(v) The pressure P is positive and finite everywhere inside the fluid 

(vi) The pressure vanishes at the fluid boundary with the vacuum. 

(vii) The speed of sound Vs = {dP/dpY^"^ is causal (0 < u < 1). 

(viii) The speed of sound = {dP/dpY^^ monotonically decreases with increasing r. 
(ix) Both the pressure and mass density are decreasing functions of r: dP/ dr < and 



This list is not necessarily exhaustive but does include all of the requirements 
presented by Delgaty and Lake in their extensive review of exact fluid sphere solutions. 
Supplementary conditions that define what Finch and Skca call "interesting solutions" 
have been reviewed in that manuscript. These extra conditions include: an explicit 
equation of state in the form P = P{p), and non- numerical solutions to the differential 
equations that governing the behaviour of some of the dependent functions. While the 
latter condition will be met in what follows, the we have not succeeded in providing 
an equation of state outside of giving both the pressure and mass density in terms of 
functions of a radial coordinate. 

In comparing the behaviour of 127 known solutions against the requirements listed 
above, Delgaty and Lake found that only 9 of those satisfied all of the conditions. 
Clearly such solutions arc rare. Many solutions have pathologies associated with 
negative pressures, and/or non causal sound speeds. As will be shown in what follows 
many of the new solutions also demonstrate such pathological behaviour. 

The outline for this manuscript is as follows: the next section introduces the 
coordinates, the line element, and the form of the Einstein-Maxwell equations to be 
solved. A strategy for solving for the unknown functions is then presented. Section 
3 then applies the method to obtain a general form of the solutions in terms of 
analytic functions. Special cases are analyzed to determine which solutions and what 
parameters lead to physically relevant behaviour by meeting all of the conditions listed 
above simultaneously. Finally, the last section concludes with some discussions. 

2. The Einstein-Maxwell Equations 

In curvature coordinates, the line element of a spherically symmetric, static spacetime 
can be written in the form: 



polar angles and (p, and i/ and A are r-dependent functions. 

The coupled Einstein Maxwell equations can be written in the following form: 



dp/dr < 0. 





, 1 



(2) 




(3) 



(4) 
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where the subscript (r) represents a derivative with respect to r and k = SwG/c'^. In 
what follows, geometric units arc used where G = c = 1. 

The charge function q{r) is obtained by integrating the charge density cr(r) over 
a spherical proper volume: 



q{r) = ^ [ a{r)r'^e^/'^dr. 
2 Jo 



In what follows, it will be assumed that the charge function is known and this will 
be used to obtain the charge density once a solution for the metric function A(r) is 
determined. 

The spherical symmetry of the problem allows only a radial component to the 

electric field which can be computed from an electrostatic potential, (j){r). This leads 
directly to the one non-zero component of the Maxwell field tensor; 



dr r' 



2" 



Equation (2) provides a single equation for A and since the the left-hand-side is a 
total derivative of the function r(l — e~^) multiplied by one introduces the mass 
aspect: 



M{r) = ^[pr'dr+\f^'^,dr 



which leads directly to an expression for A(r): 

..e- = l-^ (5) 
r 

Therefore if the mass density and the charge aspect are given as functions of r one 
can find an expression for the mass aspect and therefore for the function z{r). 

Assuming that the metric function v has been determined one can immediately 
obtain an expression for the pressure by taking a sum of the equations (2) and (3). 
Then one obtains: 

C ^ Z Z 

k(P + p) = (Vr + Xr) = -fr ' (6) 

r r r 

Finally using equations (2) - (4) one can obtain the following equations by 
introducing the function y = e"!"^: 

Ir'^zyrr + {r^Zr - 2rz)yr + (^Zr -2z + 2- y = (7) 

2r'^zyrr + [r^Zr - 2rz)yr + {brZr + 2z-2 + 4Kpr'^)y = (8) 
2r'^zyrr + [r^Zr + &rz)yr + {rzr + 2z-2- AKPr'^)y = (9) 

While the third equation (9) couples to equation (6) the first two equations are 
equivalent and lead to second order linear differential equations for the function y. 

The procedure for constructing a solution is as follows: (i) give the mass density 
and the charge aspect functions as explicit functions in the radial coordinate r; (ii) 
compute the mass aspect which leads to an explicit solution for the function z(r); {Hi) 
obtain a solution for the function y{r) using the differential equation (7) or (8); {iv) 
finally solve for the pressure, p{r) using equation (6). 

The advantages to this procedure are that one can begin with a radial dependence 
that guarantees that some of the requirements that the solution be physically 
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significant are met from the very beginning. Here it assumed that solutions for the 
conditions on the metric functions will follow if the input is physically motivated. 
Given that there will be two constants of integration associated with the solution for 
y and some free parameters associated with the choices for a and p, there should be 
enough freedom to ensure that the remaining equations arc satisfied. 

Of course the difficulty of solving the equations for y{r) depends on the choices 
that are made for p{r) and q{r). In the remainder of this paper we present a choice that 
leads to some easily solvable equations after making the appropriate set of coordinate 
transformations. We then explore these solutions to determine whether or not they 
can meet the conditions for physical acceptability. 



3. A Set of Solutions 



One must begin with defining the mass density as a function of the radial coordinate 
and we choose: 

P = Pa- Pit'^ - P2r^ (10) 

The charge function must also be given and in order to ensure that it is concentrated 
at the outer edge of the sphere it is given as a simple power function: 

q = q2r\ (11) 

With these given, the first metric function we can solve for is A. Using the substitution 
z = e~^, we have 

, 2M(r) 

z = 1 —, 

r 

where the mass function M(r) is given by 

r 

M{r) = ^ J (^t,p+^y^dr. (12) 



We integrate and obtain 

z = 1 — ar^ — br^ — cr^, 
where a, b, and c are constants: 

a = ^Kpo; (13) 
b = ^Kpi; (14) 



and 



c=]j{ql- KP2) 



To solve for the second metric function we use the Einstein equations written as 
linear second-order differential equations for y = e"/^. In particular, where subscripts 
denote derivatives with respect to the given variable, 

2r'^zyrr + [r'^Zr - 2rz) yr + (rzr -2z + 2- 4^^ y = 0. 
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A change of variable to x = gives a mass density from (|10p and a charge 
function from (jlip of, respectively, 



and 



The differential equation is now 
1 



Kpi 1 

10 7 



(8(72 - i^P^) X 



where 

2; = 1 - aa; - - cx^ . (15) 
We choose Kp2 = Sql so that p2 > and the constant c becomes 

9 1 

c = -95 = --KP2 < (16) 
and the differential equation, with another change of variable 
dx f dx 



z J y/l — ax — bx^ — cx'^ 



can now be written as the linear second-order constant coefficient differential equation 
of a simple harmonic oscillator: 

%C + dy^O. (17) 
The constant d depends on pi: 

d = —npi = (18) 
20 4 ^ ' 

Equation ([T7]) has the following three types of solutions: 

y^C\+ C2i (19) 
y = Cie^« + C2e-^« (20) 

and 

y = Ci sin (Vdi) + C2 cos (Vde) (21) 

for d = 0, c? < 0, and d > 0, respectively. 

The final function to be found is the pressure P. We have the equation 

k{P + p) = - [Vr + K) ^ -l^r - — (22) 

r r r 

which can also be written in terms of the new coordinates x and ^: 

kP = 4^1^ - 2z^ - up. (23) 
V 

The constants Ci and C2 are determined by conditions on y and the pressure. At the 
boundary of the fluid sphere, which is determined by P (xq) = 0, the metric functions 
must match the external Reissner-Nordstrom metric, given by 

e^^e-^ = l-^^-^. (24) 
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2/(6) = (25) 

A possible second condition comes from setting the boundary where p{xf)) — 0. 
In this case, equation (|22p imphes that the derivatives of v and A with respect to r 
must be equaL This can be written as 

2y^(a;o) = z^{xo) (26) 

This second condition, however, is not necessarily required. The mass density does 
not have to vanish at the boundary. In this case the boundary conditions must be 
modified. We still impose the junction condition with the external Reissner- Nordstrom 
metric, equation (|25p . however, since the mass density does not go to zero at the 
boundary, equation (|26p no longer applies. Instead we solve for the pressure, given by 
equation ([251) . goi^^g to zero. We obtain 

V^o _ I 2za;(a;o) + kp(xq) \ 



y(Co) \ 4^z(a;o) / 



and since y{£,o) — \J z{xq) from the junction condition, we have 

yco = \ (2zi:(a;o) + ^^(xo)) . (27) 
The charge density cr can be found with q and z determined. The charge density 

is 

= (28) 
= 8\/^\/x\/z. 

In addition, the total mass of the fluid sphere can be found using the relation 
m(r) = M{r) + 

2r^ 

where M(r) is given in equation (H^). We solve and find the total mass out to the 
boundary xq to be 

»7^(a;o) = '^x^q'^ + ^4^^ (29) 

while the total charge to the boundary is 

g(a;o) = ^cxl- (30) 

Yet another condition that could be used is to match the metric functions to the 
extreme Reissner- Nordstrom solution, where we equate equations ((29| and pO)) : 



2^0^ + — V cxq (31) 

We now consider the following cases for the solutions of equation (I17|) . 
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When pi — the equation for y is hnear in ^, equation We note that if the 

mass density does not go to zero at the boundary, it is possible to have a positive 
pressure solution. That is, the contribution from the mass must be greater than from 
the charge. In this case, we use the junction condition ((25|) and equation ([27|) to solve 
for the constants Ci and C2. We find 



Ci = \/ z{xq) - C2^0, 

as before, and 

^^ = 4 + 2^°- 

Figure ^ shows a solution where the mass density is always greater than the 
charge density and the pressure starts positive and monotonically goes to zero at the 
boundary. The case shown has the following choices for a, c, and xq: 

a = 1 

1 

' = -6 

Xq ^1 



The speed of sound in the fluid sphere is also of interest. This speed should be 
positive and causal, that is 

dP , , 

< — < 1. (32 
dp 

Figure ([1]) is a plot of the speed of sound in this medium which indicates that the 
speed is certainly not causal in this case. 

If this solution is matched to the extreme Reissner-Nordstrom solution, we impose 
the condition given in equation ((3T|) . Where 6 = 0, the boundary xq is solved for: 

1 a2 
2^0 = 7:- 



2(-c)- 

Unfortunately, this radial boundary also produces negative pressures over part of the 
fluid sphere. Figure ([2]) shows the example with 

a = 1 

1 

Xo = 1 

The speed of sound is singular at the center of the sphere and decreases when 
both the pressure and mass density are decreasing The speed is negative when the 
pressure increases. The speed is plotted in Figure 

If we impose the condition p{xq) = 0, the boundary of the fluid sphere is 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 08 1 

radial coordinate, x radiai coordinate, x 



Figure 1. Plots of: (a) the metric functions %p and z; (b) the pressure P, the 
mass and charge densities Kp and Ktr; (c) the mass and charge functions m(r) and 
q{r) and (d) dP/dp for fe = 0. 



and equation ([26)1 implies 
16 

The junction condition (|25p produces 

a 



Ci = y/z{xo) ~ —Co- 
in this case, the pressure, from equation (|23p. reads 

kP = 4\/i ( 7^-%7-: ) - a - 2ca;^ (33) 



The pressure will be positive at the center (x = = 0) so long as Ci < \. 
Unfortunately, having the two conditions on the pressure, that the pressure is positive 
at the center and vanishes at the boundary, is not sufficient to prove that this is 
a positive pressure case. Numerical evaluation of the pressure shows that it has an 
additional point between the center and the boundary at which it goes to zero. That is, 
the pressure starts positive, then becomes negative before vanishing at the boundary. 
We illustrate a particular example of this behavior in Figure with the following 
values for the constants a and c: 

a = 1 

1 
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0.5 1 1.5 

radial coordinate, x 




radial coordinate, x 



Figure 2. Plots of: (a) the metric functions and z; (b) the pressure P, the 
mass and charge densities Kp and Ktr; (c) the mass and charge functions m(r) and 
g(r) and (d) dP/dp for the extreme Reissner-Nordstrom solution and b = 0. 



3.2. Case d < (b > 0, pi < 0) 



This case involves the real exponential solution for y given by equation (|20p . If we 
do not impose the condition that the mass density goes to zero at the boundary, we 
can solve for the constants Ci and C2 using the junction condition (|^ and equation 
(|27p . These two equations result in 



2 



V4 



4 + 4-0 + 2^0 



Co 



These constants do give positive pressures as shown in the following example. Choose 
the following values for the constants a, b, c, and xq: 



a = i 
b = 1 
c = -2 
a;o = 0.09 
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0.5 1 1.5 

radial coordinate, x 




radiai coordinate, x radiai coordinate, x 



Figure 3. Plots of: (a) the metric functions and z; (b) tlie pressure P, the 
mass and charge densities Kp and kct; (c) the mass and charge functions m(r) and 
q{r) and (d) dP/dp for 6 = and p{xo) = 0. 



Our results are shown in Figure ([4]). 

Once again, the speed of sound in this medium ^ is evaluated and, as shown 
in Figure (|3]), is negative for this case. This is due to the fact that the mass density 
increases over the range of interest while the pressure decreases. 

A match to the extreme Reissner-Nordstrom solution leads to a boundary given 
by the solution to equation pTjl : 

So = f — ^ ± -\/-c - ab 

where the minus sign case is chosen to ensure that the mass density remains positive 
over the region of the fluid sphere. Using the same values for the other constants as 
the previous example, Figure ([5]) is obtained. 

The speed of sound becomes singular when dp = 0. 

If the mass density is forced to go to zero at the boundary of the fluid sphere 
{p{xq) = 0), we solve for the boundary: 

xo^ + -^\IpI + ^PoP2 (34) 

2/52 2p2 V 

hh 1 



16(-c) ' 16(-c)^^^^^ + ^^«(-^) 
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Figure 4. Plots of: (a) the metric functions and z; (b) the pressure P, the 
mass and charge densities Kp and kct; (c) the mass and charge functions m(r) and 
q{r) and (d) dP/dp for fe > 0. 



The two conditions given by equations 
Ci and C2, 

1 



2 



\/ z{xq) + 



and (pS)) allow us to solve for the constants 

Zx{xo) 



C2 = -e^ 
2 



V ^(2^0) 



2V-d 



The pressure, from equation ((23l) . is now 



26 



— a — bx ~ 2cx 



(35) 



yCieV-'^? + C2e 

The condition for the pressure to be positive at the center of the fluid sphere is 
C\-C2 a 
C1 + C2 2Vb' 

This clearly requires Ci > C2 but this sets no clear condition on (xq) due to the 
presence of the real exponentials in Ci and C2. The pressures appear to be negative 
when the contribution from the charge overwhelms the contribution from the mass, 
especially near the boundary where the mass density is virtually nil. The constants 
a, b, and c are the same as in the previous example while xq is now given by equation 
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(a) 




radial coordinate, x 




radiai coordinate, x 



(b) 



charge 
-emass 
pressure 



0.05 0.1 0.15 

radiai coordinate, x 



0.05 0.1 0.15 

radial coordinate, x 



(d) 



Figure 5. Plots of: (a) the metric functions and z; (b) tlie pressure P, tlie 
mass and cliarge densities Kp and kct; (c) the mass and charge functions m(r) and 
g(r) and (d) dP/dp for extreme Reissner-Nordstrom solution and fe > 0. 



(|34p . Figure ([6]) shows this negative pressure case. The speed of sound in this case 
becomes singular at the point of maximum mass density (where dp — 0). 

We can also investigate the sub-case of po = with both conditions on p{xo). 
Thus, both the mass and charge densities are zero at the center of the fluid sphere. 
This is the only relevant case in which to examine this sub-case as the mass density 
will remain positive for a large enough choice of pi. When p{xo) > we desire the 
mass density to be large relative to the charge density, so the following values for the 
constants a, 5, c, and xq are chosen: 

a = 
6 = 10 



Xq = 0.3 

We do in fact obtain positive pressures as shown in Figure The speed of sound 
in this case remains negative due to the increasing mass density with a decreasing 
pressure. The speed is shown in Figure ([7]). 

The extreme Reissner-Nordstrom solution has a boundary given from equation 
dSIl) of: 

4(-c) 



xo 



62 
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0.2 0.4 0.6 0.8 

radial coordinate, x 



Figure 6. Plots of: (a) the metric functions %p and z; (b) the pressure P, the 
mass and charge densities Kp and kct\ (c) the mass and charge functions m(r) and 
q{r) and (d) dP/dp for fe > and p{xo) = 0. 

This leads to negative pressures over part of the sphere as in the example given in 
Figure ^ where the following values for the constants are used 

a = 
c = - 2 

The speed of sound in this case begins positive yet becomes negative as the pressure 
increases and then decreases. This speed is given in Figure ([8]). 

However, when p{xo) — and po = the mass density remains below the charge 
density and, hence, negative pressures result. For example, we choose the following 
values for b and c: 

b = 2 
c = -2. 



The resulting pressure is shown in Figure 

The speed of sound in this case becomes singular when the mass density peaks 
(i.e. dp = 0). 
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Figure 7. Plots of: (a) the metric functions and z; (b) the pressure P, the 
mass and charge densities Kp and kcf; (c) the mass and charge functions m(r) and 
q{r) and (d) dP/dp for a = and b > 0. 



3.3. Case d > (b < 0, pi > 0) 



This case involves the sine and cosine solution for y given by equation (|2ip . Now, if we 
use the junction condition (|25p and equation (|?7|) . the constants Ci and C2 become 



Ci y/ z{xa) sir 



1 /a 

75 U^" 



Xq 



CXq 



Ci = v/z(a:o)cos V 7^0 - ^ + o 



\/6 V2 



sm - 



1^" 



We obtain positive pressures in this case as well. The following example, Figure((TU]), 
uses the following values for the constants a, 6, c, and a;o- 

a = 1 
h = -\ 
c = -2 
a;o = 0.05 

The speed of sound in this shown in Figure pU)) , does follow the conditions 

given by equation 
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0.05 0.1 0.15 0.2 

radial coordinate, x 



Figure 8. Plots of: (a) the metric functions and z; (b) the pressure P, the 
mass and charge densities Kp and Ktr; (c) the mass and charge functions m(r) and 
g(r) and (d) dP/dp for extreme Reissner-Nordstrom solution, a = 0, and b > 0. 



Thus, this case satisfies the conditions for a physically realistic solution, that: at 
the center of the sphere, the mass m vanishes and e'^ = 1; the metric functions 
and z match the external Reissner-Nordstrom metric at the boundary; the pressure is 
positive and monotonically goes to zero at the boundary; the mass density is positive; 
and the speed of sound in this medium is both positive and causal. 

Matching to the extreme Reissner-Nordstrom solution results in another situation 
where the pressure does not monotonically vanish at the boundary, given in Figure 
(jlll) for the same values as in the last example. The boundary is given from equation 

The speed of sound is fractional but does not remain positive as illustrated in 
Figure (fTTj) . 

The condition that the mass density goes to zero yields the boundary 

Pi 1 n> ] 

2p2 2p2 * 

1 
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Figure 9. Plots of: (a) the metric functions •Lp and z; (b) tlie pressure P, tlie 
mass and cliarge densities Kp and kct; (c) the mass and charge functions m(r) and 
q{r) and (d) dP/dp for a = 0, b > 0, and p{xo) = 0. 



The junction condition and equation (|26p are used to solve for the constants Ci 
and C2: 

C\ = \/ z{xo) sin VdCo + ^''^^°} cos Vd^o 

C2 = \/ z{xo) cos - ^'^^^'^ sin \/d^o- 

2v d 

The pressure, from equation (|23p . is given by 

yCisinVd^ + C2Cos%/de/ 
The condition that the pressure be positive at the center leads to the relation 
Ci ^ a 



C2 AVd' 

which is not satisfied, in general. We show an example in Figure ([T2|)with negative 
pressures where the following values were used 

The speed of sound in this case, begins at a positive fractional value but decreases 
through zero to negative values as the pressure first decreases and then increases to 
zero at the boundary. 

Finally One might ask which solutions (of the 127 
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radial coordinate, x radial coordinate, x 



Figure 10. Plots of: (a) the metric functions and z; (b) the pressure P, the 
mass and charge densities up and kit; (c) the mass and charge functions m(r) and 
q{r) and (d) dP/dp for 6 < 0. 



4. Conclusions 

It has been shown that by giving the mass density and charge aspect functions 
exphcitly was fourth-order polynomials in the curvature radial coordinate leads to a 
simple linear second-order differential equation with constant coefficients that can be 
easily solved. With an appropriate choice of the constants appearing in the polynomial, 
some of the solutions can be made to satisfy the requirements for begin physically 
acceptable. Given the small number of known solutions for both c;liarged and neutral 
fluid spheres that obey such criteria, the configurations of charge and mass found 
in these solutions may be expected to form from a realistic collapse to a condensed 
spherical object. 

While matching to the external Reissner-Nordstrom solution causes no problems 
since this only determines the values for the constants of integration, it has been 
shown that with this family of solutions it is impossible to find realistic configurations 
that lead to extreme Reissner-Nordstrom solutions. In all cases where the solution 
was matched to the extreme Reissner-Nordstrom vacuum solution, problems arose in 
maintaining a speed of sound that is bounded by the speed of light. 

It remains to be determined how stable these configurations are, and this will 
be studied in the future. Given that the equation of state for the fluid has not been 
written in the form p = p{p) does make the stability analysis a bit more difficult since 
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Figure 11. Plots of: (a) the metric functions and z; (b) the pressure P, the 
mass and charge densities rep and kit; (c) the mass and charge functions m{r) and 
q{r) and (d) dP/dp for extreme Reissner-Nordstrom solution and 6 < 0. 



it is not known if one can use adiabaticity arguments without directly checking them 
independently. However a general stability analysis of this and other known exact 
solutions should be interesting from the point of view of trying to understand what 
charged fluid configurations could possibly be formed from the gravitational collapse 
of material that has an excess charge density. 

Finally, one should be able to extend this method to other explicit functions for 
both the mass density and charge aspect. Higher order polynomials are expected 
to produce similar results and one can expect that these would force more of the 
charge to be distributed closer to the surface of the sphere. Further studies on what 
the most general polynomials might be and the relationships between the charge and 
mass distributions will be studied in the future. 
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Figure 12. Plots of: (a) the metric functions •Lp and 2; (b) tlie pressure P, tlie 
mass and cliarge densities Kp and kcf; (c) the mass and charge functions m(r) and 
q{r) and (d) dP/dp for b < and p{xo) = 0. 
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